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Abstract 

We calculate one- and two-point correlators of winding operators in the matrix model 
of 2D string theory compactified on a circle, recently proposed for the description of string 
t^J- ■ dynamics on the 2D black hole background. 

o\ ' 
o 

O '■ 1 Introduction 

The matrix quantum mechanics (MQM) approach to the 2D string theory has been proven 
very effective, in particular, in the case of the compactification on a circle of a radius R []lj, ]!|. 
Recently, it was shown that this approach is capable to capture such important phenomena 
as the emergence of a black hole background and condensation of winding modes (vortices on 
the world sheet). 

From the statistical-mechanical point of view such a model describes a system of planar 
^ ■ rotators (XY-model) on a two-dimensional surface with fluctuating metric. In this sense it is 
a natural generalization of its analogue on the flat 2D space studied long ago by Berezinski, 
Kosterlitz and Thouless (BKT). The role of the fluctuating surfaces is played by the planar 
graphs of the MQM and the vortices appear due to the compactification of the time direction. 

As it is usual for the matrix model approach, it enables us to apply the powerful framework 
of classical integrability to the computations of physical quantities, such as partition functions 
of various genera or correlators of physical operators. 

The appropriate MQM model for the compactified 2D string is the inverted matrix oscillator 
with twisted boundary conditions 0. The partition function of such a system as a function of 
couplings of the potential for the twist angles appears to be a r-function of the Toda lattice 
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hierarchy M, [|. As such, it satisfies an infinite set of partial difference-differential Hirota 
equations. The couplings play the role of "times" of commuting flows in the Toda hierarchy. As 
was shown in || they are at the same time the couplings of vortex operators of various vorticity 
charges. The tree approximation for the string theory corresponds to the dispersionless limit of 
Toda equations. 

In this paper we will use Hirota equations to calculate all one-point and two-point correlators 
of vortex operators with arbitrary vorticity charges for the compactified 2D string in the tree-like 
approximation in the regime of the radius of compactification \Rkt < R < Rrt (Rrt is the 
BKT radius), when condensation of vortices of the lowest charge n = ±1 takes place. Unlike the 
standard BKT model in flat 2D space, where the Debay screening in the vortex plasma destroys 
the long range correlations, the BKT vortices in the fluctuating 2D metric can be described as 
a conformal matter with central charge c = 1 coupled to the metric. 

An important particular case of this model, for R = ^Rrt, claimed in || to describe the 2D 
string in the dilatonic black hole background, due to the conjecture of |5|] based on the results 
of 0, |fl. This conjecture establishes a kind of a weak-strong duality between two CFT's: the 
dilatonic black hole and the so called Sine-Liouville theory. We will try to compare our results 
with a few calculations of correlators in the continuum CFT performed in J/J |6], H and discuss 
the difficulties of the direct comparison. 

In the next section we briefly describe the basic ideas of the MQM approach to the 2D string 
theory with vortex excitations; the Toda integrable structure of the model allows to formulate 
the Hirota equations describing the dynamics of vortices. In section 3, using Hirota equations 
the two point correlators of arbitrary vorticities are expressed through the one point correlators. 
In section 4 the one point correlators are explicitly calculated. The explicit expressions for the 
two-point correlators can be found in section 5. Section 6 is devoted to the comparison with 
the results of continuous (Liouville) approach and to the conclusions. 

2 Toda hierarchy for the compactified 2D string 

2.1 Matrix model of 2D string and coupling to the windings 

In this section we recall the main results of the paper || leading to the Toda hierarchy description 
for the compactified 2D string theory. 

This theory is represented by the Polyakov string action: 

S(x, g) = — [ d 2 aJ detg[g ab d a xd b x + //], (1) 

47T J 

where g a b is a two dimensional fluctuating world sheet metric and the only bosonic field x(a) is 
compactified on a circle of a radius R: 

x{a) ~ x(a) + 2tcR. (2) 

Such a field admits in general Berezinski-Kosterlitz-Thouless vortices (windings) in its configu- 
rations. If we denote by x n±1,n±2 '"'(a) a configuration of the bosonic field containing n±i vortices 
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of vorticity charge ±1, n± 2 vortices of vorticity charge ±2, etc, the partition function of the 
theory in this sector will be 

Z n±un±2 ,,. = I Dg(a) I !);■" ■ ■■■ -■(o)< (3) 



Obviously, only the total vorticity YlkL-oc ^ n k is conserved by topological reasons, but we will 
distinguish configurations with different individual vortices as point-like objects with vorticities 
k, to give a statistical-mechanical meaning to the system. 

Instead of fixing the numbers of vortices we can also introduce the partition function with 
fixed fugacities £±i>£±2> • • • of the vortices with given charges: 

Zt±i,t±2,- — X! ^fci^if 2 ' ' ' Zn±i,n± 2 ,-- (4) 

n±l,™±2,--- 

It is well known that in the case of flat 2D space (g a b = 5 a b) such a system can be described 
in a dual, coulomb gas picture, in terms of a 2D field theory of one massless scalar field x(a) 
(not necessarily compactified) perturbed by vortex operators built out of this field. In the case 
of fluctuating metric this scalar field gets coupled to the Liouville field. The action (at least in 
the small coupling regime <fi — > oo) looks as 

S = -L [d 2 z [(dx) 2 + (d(j)) 2 - 47^0 + //0e- 2 * + V x'e^ R - 2 ^e inR£ ), (5) 

where x = xl + xr, x = xl — xr (xl and xr are left and right movers). The total curvature 
is normalized as / d 2 z ^fgR, = 2n(2 — 2h), where h is the genus of the surface. Note that we 
chose the central charge c = 26 and the Liouville "charges" of vortex operators to make them 
marginal. The compactification radius R remains a free parameter of the theory. If we choose 
the couplings as 

A'^^i+M. (6) 
the model will turn into the Sine-Gordon theory coupled to 2D gravity: 

S SGG = -L fd 2 z [(Ox) 2 + (d<f)) 2 - 47^> + /i'0e- 2 ^ + X'e iR - 2 ^ cos Rx] . (7) 
Ait J 

As its counterparts on the flat space, this theory is supposed to be in the same universality class 
as the compactified 2D string theory ([!]). 

In particular case of vanishing cosmological constant fi' and R = 3/2 (here Rkt = 2) the 
CFT (0) was conjectured to be dual to the coset CFT describing the string theory on the 
dilatonic black hole background [|, §] (see for the details). 

The matrix model version of this theory is represented by the partition function of the MQM 
on a "time" circle of the radius R 

Zn[\ g]= I W} SU (N) exp I £ A n trfi n Z N (Q, g), (8) 
J Vnez / 
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where 

Z N (n,g)= I VM(x)e~ tr C R H^M^ + v(M)}_ (g) 

JM(27rR)=ntM(0)O 

We take the twisted boundary conditions for the hermitian matrix field M(x), where the twist 
matrix Q can be chosen without lack of generality as a Cartan element of the SU(N) group 
Q = diag(zi, z 2 , ■ ■ ■ , z N ). The matrix potential is, for example, V(M) = \M 2 — ^-^M 3 . The 
couplings A ra will play the same role here as A^ in (|5|) although the exact relation between them 
depends on the regularization procedure (the shape of the potential V(M)). 

Following the usual logic of the "old" matrix models, the Feynman expansion of this model 
can be shown (see || for the details) to describe the compactified lattice scalar field X{ living in 
the vertices i = 1, 2, 3, . . . of 3 type planar graphs. This field can have vortex configurations: 
the vortices of charges n = ±1, ±2, . . . occur on the faces of planar graphs and are weighted 
with the factors A n . So the model (H) represents a lattice analogue of the continuous partition 
function for the models (|I]), (|4|) or The sum over planar graphs represents the functional 
integral over two dimensional metrics of the world sheet and due to the standard 't Hooft 
argument 1/N expansion goes over the powers 1/N 2 ~ 2h where h is a genus of the world sheet. 

This model has a long history. In [[Uj some estimates of the matrix approach were compared 
to the predictions of (0) and appeared to be in full agreement. In || the model (P was 
formulated in terms of a twisted inverted matrix oscillator (||) . The Toda hierarchy description 
of the model (for the usual stable oscillatorial potential) was proposed in [||. In the paper 
TT] the dual version of the matrix model of [|I|, [RJ, the Sine-Gordon model on random graphs, 



was used to find (or rather correctly conjecture) the free energy of (|5]) as a function of the 
cosmological constant /i', Sine-Gordon coupling A' and radius R. 

The twisted inverted matrix oscillator of appeared to be very useful for the identification 
of the free energy of the whole theory (^) and of the corresponding matrix model (^) with the 
r-function of the Toda chain hierarchy of integrable PDF's ||. To see how this description 
emerges we pass to the grand canonical partition function: 

oo 

Z,[\,g}=Y,e- 2nRt " N Z N [X,g]. (10) 

N=0 

The chemical potential /i happens to be the cosmological constant similar to /2 appearing in (|7|) 
(see |]12[ and for the details). Then following the double scaling prescription for the c = 1 



matrix model |13j corresponding to the continuous limit for the lattice world sheets we should 



send in (H) N — > oo and g g cr u in such a way that at the saddle point of the sum over N in 

» = 2^RM l0gZN[X > 9] 
remains fixed. This amounts to shifting M — > M + ^j- and neglecting the cubic term in the 
potential V(M), or rather treating it as a cut-off wall at M — * oo, which gives V(M) ~ — \M 2 . 
Then the integral |2j can be calculated exactly, giving: 

Z-j — Zai _ . 

= m / n ^— {ql/2 _ q - 1/2) n q i/2 Zj _ q -i/2 Zjl - en) 
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where u(z) = | J2n X n z n , q = e 2mR . 

For the identification with r-function it is convenient to redefine the vortex couplings as 



t = — (12) 

ln q n/2_ g -n/2- ^> 



Then plugging (|Tl"D into (|10D we recognize in it the r-function of Toda lattice hierarchy JT^, |15 



n [t] = e"E» (9* e^T-ZivM = e-Sn^-^M, (13) 

Af=0 

where the charge / giving an extra "lattice" dimension to the Toda equations appears to be an 
imaginary integer shift of the cosmological constant ix. The couplings t n of vortices, together 
with to = fx turn out to be the "times" of commuting flows of the hierarchy. Due to this the whole 
sector of the theory describing the dynamics of winding modes of the model ([!]) is completely 
solvable and the calculations of particular physical quantities, such as the correlators of vortex 
operators, are greatly simplified. In particular, the r function (|T3| ) represents the generating 
function for such correlators in the theory fl7|): 



gn 

dX^ ■ ■ ■ d\i 



fch-in = Er^ lo g T 



(14) 

A±2=A±3=---=0 



with [x and A±i fixed. 

Using the Hirota equations of Toda hierarchy we will calculate in this paper the one-point and 
two-point functions of vortex operators of arbitrary vorticities in the spherical approximation. 

2.2 Hirota equations 



One can show |Tb| that the ensemble of the r-functions of the Toda hierarchy with different 
charges satisfies a set of bilinear partial differential equations known as Hirota equations which 
are written as follows (the derivatives are taken with respect to t n rather than A n ) 

§ Pj+i(-2y+)pj{D + ) exp f £ y k D k I r i+l+1 [t] ■ n[t] = 
j=0 \k^o J 

§ Pj-i{- 2 y-)Pj(D-) exp f £ y k D k I r i+l [t] ■ r m [t], (15) 
3=0 Wo / 



where 



y± = (y±i,y±2,y±3, ■ ■ •)> ( 16 ) 

At = (D ±1 ,D ±2 /2,D ±3 /3,...) (17) 



represent the Hirota's bilinear operators 



D n f[t] ■ g[t\ = ^-f(t n + x)g{t n - x) 



x=0 
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and pj are Schur polynomials defined by 

oo / oo \ 

J2Pk[t]x k = ex V ]Tt n x n . (19) 

fc=0 \fe=l / 

Since the r-function of the Toda hierarchy is related to the free energy by 

T a \ji,t]=exp(F(jM-is)), (20) 
the Hirota equations lead to a triangular system of nonlinear difference-differential equations 



for the free energy of the model (|10T)-(|lll). Due to CpL4[) it can be actually thought of as a system 



of equations for the correlators which we are interested in. 
2.3 Scaling of winding operators 

It turns out that there is a scale in the model which corresponds to the scale given by the 
string coupling constant in the string theory. Since g s ~ e < ^ > it can be associated with the 
Liouville field or the cosmo logical constant fi' coupled with it. All couplings X' n have definite 
dimensions with respect to this scale. In the Toda hierarchy the counterparts of fj,' and X' n are 
\i and t n correspondingly. As one can immediately see from (|5|) by shifting the zero mode of the 
Liouville field, the corresponding dimensions of the couplings t n with respect to the rescaling of 
the cosmological constant \i are 

A[*„] = 1 - (21) 

This scaling leads to an expansion for the free energy which can be interpreted as its genus 
expansion 

oo 

F = H F h, (22) 

h=0 

F = r 2 /o(^) + !/i 2 log£, (23) 
F h = e h - 2 fh(uj;s), h>0, (24) 
where the quantities fh are the functions of the dimensionless parameters 

S = {8 2 , S 3 , • • •), S n = I (-^)" /2 e [tn] tn, (25) 

u = fi$, with f = (X^R^iy^R = {txt^R- 1))^. (26) 

Remarkably, the Toda equations (|T5| ) are compatible with this scaling. 

The spherical limit of the string partition function is described by the dispersionless Toda 
hierarchy which can be obtained taking £ — ► 0. In this limit the explicit expression for the 



partition function with vanishing s n has been conjectured in |T1] and proven in . The result 
formulated in terms of X (u) = d^f (uj) reads 



Tt x ° - e -V x o (27) 



_ K-l 

w = e 
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or, integrating X twice over u 
F =~// 2 (Rlog£ + X )+£ 



\4R-1 + 4 R-l 



-Xo 



(28) 



This result will serve as an input for the following calculation of the correlators in the spherical 
limit. 



3 Two-point correlators from the Toda hierarchy 

3.1 Equation for the generating function of two-point correlators 



To find equations for the two-point correlators in the spherical limit let us take in ( T5|) i — 
and the coefficient in front of y n y m , n,m > 0. Then we obtain the following equation 

(D + ) — 2p n (D + )D m — 2p m (D + )D n + D n D m \ r s+1 ■ r s = D n D m r s ■ t s+1 , 

2p n+m (D + ) - p n (D + )D m - p m (D + )Dn\ t s+1 ■ r s = 0. (29) 

This equation is valid also at m = or n = if we set Dq = 1. Such an equation can be 
obtained from the coefficient in front of y n in (|15|) . 

If we multiply the equation (E0) by x n y m and sum over all n and m we obtain 



2 E x n y m p n+m (D+)- E y m D m exp\j: x n D n 



- E x n D n exp E y m D m r s+1 ■ r s = 0. 

n=0 Lm=l J . 

Furthermore, it is easy to see that the first term can be rewritten as follows: 

2 



(30) 



x-y 



E x k o k E y h Dk 
xe k=1 — ye k=1 



(31) 



It is easy to check the following remarkable formula 

f(D)e F ■ e G = e F e G f (D + (OF ■ 1 - 1 ■ dO)) , (32) 

which allows to express r-functions in the equation fl30|), (|3TD through the derivatives of the free 
energy: 



2 

x-y 



E x k (D k +X k ) E 2/ fe 0Dfc+X fe ) 



V 



oo \ Y,x k {b k +X k ) / oo \ El/ fc (5fc+X fc ) 

1 + E l/ m (^ m + X m ) e*=^ + 1 + E x"(D n + X n )) e*=i , (33) 

\ m=l / \ n=l / 
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where we have introduced the notation X k = ^X k and [] 

X k = d k Fi_ • 1 - 1 • d k F_i_, k^O. (34) 



The l's in the right hand side of fl33|) arise from the definition of D . 

We show in Appendix A that in the dispersionless limit only the second derivatives of the 
free energy survive. This means that we can put the second and higher derivatives of X k to 
zero: D n X k — 0, n > 1. Besides, using fl2D) we obtain in this limit 



A '» = - , R8& F « : =*^ (35) 

~ ~ Id 2 

D n X m = 2- — - — - F := X n m . (36) 

\n\ \m\ ot n ot m 

Thus the fields X nm correspond (up to a numerical factor) to the two-point correlators, whereas 
X 0tn give the one-point ones. We also define the generating function of all X's 

oo 

F(x,y)= £ x n y m X n ,m, (37) 

n,m=0 

where we have taken Xo 5 o = 0. 

With these definitions, trivial manipulations which can be found in Appendix B, lead to the 
following equation for the generating function 



x + y 
x-y 



(e^CMO _ g^Gw)) = x d x F(x, y)e^ F{v ' y) + yd y F(x, y)e^ x > x \ (38) 



3.2 Solution in terms of one-point correlators 

It turns out that the differential equation (|3~8|) is solvable. First of all, change the variables to 







X 


y 


c ^F(x,x)d% 




X 


y 



drj 

d £ = e ^]v)r_ e ^'(M)^, (39) 



It is easy to see that 



where 



dz i 



I dZ e i F(z , z) = bg x + ^ _ 
J z v 



dz 



h(x) = J — (e^ F ^ - l) (41) 



1 We can shift the index s of the free energy F s = logr s since due to ( p0| ) the resulting equation depends only 
on the difference of the indices of two r-functions. 
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is a regular function at x = 0. 

In terms of the new variables the equation reads 



d v F(x,y) 



1 x + y 



~+~ V ( -k 
e 2 

— II v 



Thus 



2 a; — ?/ 

F(x,y) = log- 



-'•'/ 



+ 9(0, 



(42) 



(43) 



(x -y) 2 

where the function g(() is fixed by the requirement for F(x, y) to be analitycal in both argu- 
ments. 

The analitycity condition together with fl4*0|) leads to the following requirements on g((): 
i) <?(C) ~ iC; ii) #(0 ~ l°gC 2 - There is only one function satisfying both of them, which is 
given by 



(/(C) = log (e< + e~< - 2) = log (4sh 2 ^ . 



(44) 



Taking together (43), ( f44"|) and ( [l"0|) we obtain the solution for the generating function of the 
two-point correlators in terms of a still unknown regular function h(x) 



F(x,y) = log 



Axy 



x s 



(x - y) 



r sh \-(h{x) — h(y) +log — 



(45) 



From ([45]) it is easy to see that /i(x) is nothing else than the generating function for the one-point 
correlators since h(x) = F(x, 0). 



3.3 Two-point correlators with vorticities of opposite signs 

Up to now we considered the two-point correlators for vorticities of positive sign only. Now we 
generalize the results for other cases. 

First of all, the case of both negative signs can be directly obtained from the previous one. 
The corresponding equation which is derived from the coefficient in front of y_ n y_ m in the Hirota 
equation looks as (EDI) 



(D_)-p n (D_)D 

—m Pm 



T~s ■ T s+ 1 



0. 



(46) 



The only difference is that t s and t s+ \ are exchanged. Due to this there is an additional minus 
in the definition of Xo 5 _ n . So we should distinguish two types of correlators 



2 MM a&t m F * n,m^0. 



(47) 
(48) 



With this definition, the generating functions of the correlators for positive and negative vor- 
ticities are 

oo 

F±(x,y)= £ x n y m X± n<±m , (49) 

n,m=0 
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Despite of this minus all equations for the case of negative vorticities being written in terms of 
X~ m are the same as for the case of positive ones. As a result F~(x,y) is given by the same 
solution ( f4"5|) as F + (x,y). The only difference between them is the generating function of the 
one-point correlators h which should be replaced in ( 45|) by h correspondingly 

To find the correlators with vorticities of opposite signs, take the coefficient in the Hirota 
equation in front of y n y-m and set i — 1. Then we obtain the following equation 



-2p n (D + )D- m T s+2 ■ r s = -2p m (DJ)D n + D- 1 D- m D n 
p n (D + )D„ m T s+2 ■ r s = p m (D-)D n T s+1 ■ r s+1 . 



T S+1 ' T S+1; 



(50) 



The last term in the first line vanishes since it always gives rise to Dk(F s+ i ■ 1 — 1 • F s+ i) = 0. In 
general, the one-point correlators appearing from r s+ k-r s are always supplied with the coefficient 
k what follows from fl20]). Due to this in the right hand side there are no one-point correlators 
whereas in the left hand side they enter with the coefficient 2. This equation is valid also for 
m = if we take Dq = —2 instead of the previous choice. However, for n = it is never fulfilled. 

The arguments similar to those of Section 1 lead to the following equation for the generating 
function 



.4 



yd y (G(x,y)-2h-(y)) 



\F+{x,x)+h+{x) 



-d x G(x,y)e 



lF-(y,y)-h-(y) 



where we have introduced 



A = exp (-^F ) = £- R e - 9 *f°, 



G(x,y) - 

Using fl4lD we can rewrite (|5l|) as follows 



E x n y m X n ,_ 

n,m=l 



(51) 

(52) 
(53) 



Ad T (xe 



h+(x) 



) d x + 



1 



d y G(x, y) = -lye 



h~{y) 



(54) 



It is clear that the solution of this equation looks as 

G{x,y) = 2 log {ye h ' {y) ) + f{-e~ h ' {y) - Axe h+{x) ). 



(55) 



The unknown function / is determined by the analitycity of G(x,y). Requiring cancelation 
of the logarithmic singularity at y = we obtain that f(z) = 2 log z, As a result we find the 
generating function 



G(x, y) = 2 log (l — Axye 



h+(x)+h-(y) 



(56) 



Let us note that the equations (59) and (BB) resemble the equations for two-point correlators 



of the two-matrix model found in 17 
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4 One-point correlators 



4.1 Equation for the generating function 

Now we should work out an equation for the one-point correlators. We can always put t\ = — t_i 
(corresponding to Ai = A_i, see ([12"D). From the equations written below for positive vorticities 
it is clear that in this case h ± (x) coincide with each other. For this reason we omit in the 
following the inessential sign label in h(x). The case of general t's can be obtained from the 
previous one by the substitution h^{x) = h {(—t-i/ti) ±l l 2 x\. It can be easily seen from (|25|) 



and is due to the total vorticity conservation. Accordingly, 

F ± (x,y) = F ((-t_i/ti) ±1/2 x, (— t-i/ti) ±1/2 2/) , 
G(x,y) = G tl= _ i _ 1 



-t-i/h 



a/2 



-t-i/h 



-1/2 



(57) 
(58) 



The equation we are looking for arises from two facts: i) we know the dependence of the free 
energy on t±\ and fi ( p7j) , (p8|) ; ii) the system of equations for the correlators is triangular. Due 
to this we can express X±i >n through Xq^ by a linear integral-differential operator. Namely, in 
terms of the generating functions we can write 



d y F+(x,0) = K {+) h(x)+Xl 
d y G(x,0) = K^hix), 



0: 



(59) 
(60) 



where the operators are to be found. The last term in the first line is explicitly added since 
in our notations -F + (0, 0) = Xq = 0. Moreover, we can rewrite these equations in terms of the 
generating function for the one-point correlators only. Indeed, from and (|56|) we obtain 



k i+) h{x) 

K^hix) 



1-e 

2Axe h(x) 



-h(x) 



2X, 



0,1' 



(61) 
(62) 



To find the operators we compare the derivatives of the free energy with respect to fi 
and t±i. From (23), (p5 ) and (|26|) we obtain for \n\ > 1 



d_ (d_ 

dfi \dt n 



it 



hr d 
i 







du \ ds, 



/o(w; s) 











dt±i \dt r 



(\n\+l)R 

€ 2 



These relations together with 



R 

) and 



uj- 



d_ 



1 



5fc=0 

R\n\ 



± 



{R-2)n 



' d_ 

OSr, 



UJ 



2 2 
imply that the operators are given by 
d w 



fo(y; s ) 



0J 



(l + x 



dx' 



du 



(63) 

(64) 

=o 



(65) 
(66) 



11 



where a 



2-R 



£ 2 . This expression is also valid for the case n = ±1 where the derivatives 



of the free energy contain additional terms appearing from the logarithmic term in (|23f) . 

Inserting ( |65] ) and fl6"6|) into (^1|) and fl62| ) and taking the derivative with respect to uj we 
obtain two differential equations for the generating function of the one-point correlators 



2 e -fc(x;w)_?_ 

9a; 



\ OX OUJ OUJ 



h(x; uj) 
h(x; uj) 



ouj 



(67) 
(68) 



4.2 Solution 



Let us take advantage of having two equations for one quantity. Using (|6"7f) and (|68D we can 
exclude the derivative with respect to x. The resulting equation is 

(a(2 - R)uj + 2 e — - 2(i2 - \)AxS\ d w h = 2d u Xf tl + 2(R - l)d w Axe h . (69) 

From (|7|), ([23|) and (|28]) one can obtain 



a 2 



2-R 



i? 



(70) 

Noting that the equation ( E5p is homogeneous in the derivative <9^ we can change it by <9x - 
Then after the substitution of the definitions of all entries the common multiplier 



2y/R=lC R/2 (e-^ Xo + VR^lC R/2 xe h - x « 



can be canceled. As a result we obtain the following simple equation 

(l 



^ 2 -e^ Xo e- h )d Xo h = l. 



(71) 



Xy/R^T 

This equation can be trivially integrated giving after the proper choice of the integration constant 

(72) 



e * h - ze h = 1, 



t-R/2 R-l 

where z = ■ k 



m. 



xe 



r x °. As it is easy to check this solution satisfies both equations (pTj) and 



4.3 One-point correlators of fixed vorticities 

Explicit expressions for the one-point correlators are given by the coefficients of the expansion 
in x by the generating function h(x) satisfying equation ([721). To find them we can use equation 
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5.2.13.30 in 

1 6 ^ z n T(na + b) 

b ^ n\T(n(a-l) + b+lY 1 1 

s = l + zs a . (74) 

Taking s = e^ h , a = R we obtain for b = kR 

e kh -kRY *" T « n + W (75) 

" h^.r((n + k)R-n + iy 

Also the limit b — > gives 

This means that the derivative with respect to /i of the one-point correlator is 

d 2 „ zlYnE + l) r n/?/2 „«-iy 

F = -n^-X ^7) 

dfidt n n\T(n(R-l) + l)(R-l) n / 2 ' 1 ; 

After integration over fi and taking into account the relation (|l^) we obtain the following one- 
point correlators of operators of vorticity n in the spherical limit 

/ ^ -4 \ nR+2 / n(fl-l) + l v /„ , -, \ -R-l v \ 

1 Y{nR + 1) i — — (e « — x ° _ e - {n+1) — x ° \ 

n ~ 2 sin 7mRn\T(n(R- 1) + 1)(R- l) n l 2 \ n(R-\) + 1 n + 1 J' ^ ' 

(Here the lower limit of the integration is chosen to be X = +oo. Only this choice reproduces 
B^Fq from (|28|) .) In particular, in the limit \i — > (X — > 0) corresponding to string theory 
on the black hole background (in fact, only the point R = 3/2 is supposed to describe the 
conventional black hole) we find 

K 1 (2-R)nT(nR+l) 

n 2smirnR(R-l) n / 2 (n + l)\T(n(R-l) + 2y K ' 

5 Results for two-point correlators 

5.1 Correlators with vorticities of opposite signs 

These correlators are given by the expansion coefficients of the generating function (|5*B|). This 
function can be easily expanded in e h for which the series fl75j) can be used. This leads to the 
following expansion 

oo r>k oo oo 

D(x,y) = -2j2 1 ~z k (x)z k (y) £ £ C k C k m z n {x)z m {y) 

k=l n=0 m=0 

oo min(n,m) r>fc 

= -2E z n {*)z m {y) E T C k n _ k C k m _ k , (80) 

n,m=l k=l 
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where 



B = A(R - 1)£V^*° = (R- ljeT 1 ', (81 



C k kR T((n + k)R) 

n\ T(n(R-l) + kR+l)' { ' 

The two-point correlators appear to be 

T(nR + l)T(mR + l) t 1 ^ -(n +m )^i 
2sin7mi2 2sin7rmi2 (i? - 1)(«)/ 2 

min ^ m) fc(g - i) 



fc=i 



(n - fc)!(m - k)\T(n(R - 1) + k + l)r(m(i? - 1) + jfe + 1) ' 



To get the result in the black hole limit it is enough to set Xq = in the above expression. 
It is worth to give also the result for the ratios of correlators which do not depend on possible 



leg-factors — wave function renormalisations of the operators (see for ex. [12 1): 



/C„,_ m e mil g m) k(R - l) k (n + l)!(m + l)W(n(R - 1) + 2)T(m(R - 1) + 2) 



/C n /C_ m (2-R) 2 nm ^ (n - k)\{m - k)W(n(R - 1) + jfe + l)T(m(R - 1) + k + 1) ' 

(84) 

Since they can in principle be compared with CFT calculations in the black hole limit we restrict 
ourselves only to this particular case. 

5.2 Correlators with vorticities of same signs 

To investigate this case we should expand (fi5| ) in x and y. The expansion in the first argument 
gives for n > 

9 ( 1 n 2 \ 

l K ~~ X J m=0 k=l K 

r, /-i n n i \ 

+ ^ --E^)E^ ■ (85) 

y \ n m =l fe =m K / 

The last term in the right hand side is singular and we know from the analitycity of F(x,y) 
that it should vanish. The remaining part gives the values of the two-point correlators 



T(nR + l)T(mR + l) ^ + m )R /2 
2waimR 2 shiTrmi? (R - l)(™+™)/2 
n k 
X ^ (n - k)\{m + k)\T{n{R - 1) + k + l)T(m(R - 1) - k + 1) ' ( ^°' 
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As in the previous case the black hole limit of ([SB]) is obtained vanishing X Q . The normalized 
correlators are 



JCn, n e A k(n + l)\(m + l)\V(n(R-l) + 2)V(m(R-l) + 2) 



/C„/C m (2 - i?) 2 nm ^ (n - fc)!(m + fe)!r(n(i2 - 1) + k + l)T(m(R - 1) - k + 1) 

It is clear from the definition that (86) as well as (jSTD should be symmetric in m and n. 
However, our results do not possess this symmetry explicitly. The reason is that the expansion 
of (f45|) has been done in a nonsymmetric way. It may be possible to symmetrize them by means 
of identities between T-functions like 

n j 

E r^LkCk-m = 0, 1 < m < n. (88) 

k=m 

This identity follows from vanishing of the last term in (B5|). 



6 Discussion 



The main result of this paper is the calculation of one- and two-point vorticity (winding) cor- 
relators given by the generating functions (|45D , (|56"D and (|72|) . What kind of physics do they 
describe? 

We can look at our model from two points of view: as a statistical-mechanical system 
describing a gas of vortices and clS db string theory. 



Statistical mechanical picture 

In the statistical mechanical picture we interpret the planar graphs of the MQM (|||) as random 
dynamical lattices populated by Berezinski-Kosterlitz-Thouless (BKZ) vortices. 

Our one point correlator K n ([78| ) should be in principle proportional to the probability to 
find in the system a vortex of a given vorticity n. The probability should be a positive quantity. 
However, we see that the explicit expressions ([78|)-(|79|) contain the sign-changing factors sin l Rn ■ 
Their origin is not completely clear to us but we see two possible explanations for them: 

i) They appeared as a result of wave function renormalisation in the continuous (double 
scaling) limit and should be absorbed into the "leg-factors" . The rest of the one-point correlator 
is strictly positive in the interval 1 < R < 2 which we pretend to describe. This point of view 
is supported by the fact that these factors appear as the same wave-function renormalisations 
also in the two point correlators (|83|) ,(|86|). Since their origin is completely due to the change of 
variables (fP2"D these factors appear to be a general feature for all multipoint correlators. 

ii) The sign-changing might be a consequence of the fact that our approach based on the 
Hirota equations (describing directly the double scaling limit) does not keep track of non- 
universal terms due to the ultraviolet (lattice) cut-off in the system. Such UV divergences are 
well known in the Coulomb gas description of the BKT system: they usually correspond to 
the proper energy of a vortex. In that case the correlators we calculated above are only the 
universal parts of the full correlators containing also big positive non-universal contributions. 
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These universal parts cannot be interpreted as full probabilities and have no reason to be 
positively defined. 

To prove or reject these explanations we have to see what happens with the individual 
vortices in the process of taking the double scaling limit, from planar graphs picture to the 
inverted oscillator description. It is not an easy, although not hopeless, task which we leave for 
the future. 



String theory picture 

The MQM describes also the two-dimensional bosonic string theory in compact imaginary time 
(i.e. at finite temperature). Due to the FZZ H conjecture, at the vanishing cosmological 
constant it can be also interpreted as a string theory in the Euclidean 2D black hole background 
(at least at R = 3/2) M. Our one-point correlators should contain information about the 
amplitudes of emission of winding modes by the black hole, whereas the two-point correlators 
describe the S-matrix of scattering of the winding states from the tip of the cigar (or from 
the Sine-Liouville wall). However, the exact correspondence is lacking due to the absence of 
the proper normalization of the operators. In || and || some two- and three-point correlators 
were computed in the CFT approach to this theory. It would be interesting to compare their 
results with the correlators calculated in this paper from the MQM approach. However there 
are immediate obstacles to this comparison. 

First of all, these authors do not give any results for the one point functions of vortices. 
In the conformal theory such functions are normally zero since the vortex operators have the 
dimension one. But in the string theory we calculate the averages of a type < / d 2 zV n (z) > 
integrated over the parameterization space. They are already quantities of zero dimension, 
and the formal integration leads to the ambiguity * oo which should in general give a finite 
result. Another possible reason for vanishing of the one-point correlator could be the additional 
infinite W-symmetry found in [7| for the CFT (0) (at R = 3/2, \il = 0) Q. The generators 
of this symmetry do not commute with the vortex operators V n , n / ±1. Hence its vacuum 
average should be zero, unless there is a singlet component under this symmetry in it. Our 
results suggest the presence of such a singlet component. 

Note that the one-point functions were calculated in a non-conformal field theory with 
the Lagrangian: 



4:71 



{dxf + (<90) 2 - 4R(f) + m (e~3* + e^) 



( 3 \ 
cosl-x 

K 2 ' 



*9) 



This theory coincides with the CFT (0) (again at R = 3/2 and // = 0) in the limit m — > 0, 
0o - ► — oo, with A = me"2*° fixed, where 0o is a shift of the zero mode of the Liouville field 
4>{z). So we can try to perform this limit in the calculated one-point functions directly. As a 
result we obtain JC^ ~ m 2 A 3ri_4 . The coefficient we omitted is given by a complicated integral 
which we cannot perform explicitly. It is important that it does not depend on the couplings 
and is purely numerical. We see that, remarkably, the vanishing mass parameter enters in a 
constant power which is tempting to associate with the measure d 2 z of integration. Moreover, 



2 We thank A.Zamolodchikov for this comment 
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the scaling in A is the same as in (|79|) (at R = 3/2) up to the n-independent factor A -8 . All 
these n-independent factors disappear if we consider the correlators normalized with respect to 
/C^ which is definitely nonzero. They behave like ~ A 3 ^ 1 ) what coincides with the MQM 
result. 

The two point correlators for the CFT (|7|) with // = are calculated in |J only in the case of 
vortex operators of opposite and equal by modulo vorticities. They should be in principle com- 
pared to our /C n _ n correlators from fl83"|). However, we should compare only properly normalized 
quantities, knowing that the matrix model correlators generally differ from the CFT correlators 
by the leg-factors. To fix the normalization we have to compare the quantities similar to (|S3|) 
and (|87|) which are not available in the CFT approach. 

A lot of interesting physical information is contained in multipoint correlators (3-point and 
more) which we did not consider in this paper and some of which are calculated in the CFT ap- 
proach of || |19| . The direct calculations from the Hirota equations look very tedious. Hopefully 
the collective field theory approach of [^UJ can help on this way. 

A lot is to be done to learn how to extract information about the black hole physics from the 
correlators. First of all, we have only calculated some vortex (winding modes) correlators in our 
MQM approach. The vertex (momentum modes, or tachyon) correlators are unavailable in the 
To da hierarchy approach used here and do not seem to fit any known integrable structure in the 
MQM, although they are in principle calculable directly from the MQM in the large N limit. 
On the other hand, they bear even a more important information about the system than the 
winding correlators: they describe the S-matrix of scattering of tachyons off the black hole and 
can be used to "see" the black hole background explicitly. The one-point tachyon correlators 
may give information about the black hole radiation. In general, the MQM approach seems to 
be a good chance to work out a truly microscopic picture for the 2D black hole physics. 
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Appendix A 



Let us investigate which terms will survive in the dispersionless limit in the equation fl29"|) . This 
limit corresponds to £ — ► 0. It is clear that all terms have the form 



Li=i 



T s +l ■ T s 



(90) 
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with n = kj fixed. When we rewrite them in terms of the free energy F s = logr s we obtain 

3=1 

a set of terms of a kind 



-,F S +F S . 



n 

3=1 



1j Q 

i=i j 



where E qj = I. In the dispersionless limit ( 
3=1 

F s+ i + F s 



(F s+1 + 1 , 

and (120) imply 



(91) 



2F n 



<^-F , 

OOJ 



(92) 
(93) 



The main contribution in the dispersionless limit comes from the terms with the minimal degree 
of £. The total degree of £ for a given term is 



A [t kj ] -2p + r = l-2p + r-—, 

3=1 2 



(94) 



where r is the number of odd qj. It is easy to see that the minimum is achieved when all qj = 1 
or 2 and it equals — ^ what is independent on all particular parameters and hence is similar 
for all terms in the equation. It means that only two- and one-point correlators survive. 



Appendix B 

Combine together all terms with Dk and separately with Xk in the exponents of the eq. (|33|). 
Then due to vanishing of D 2 Xk we can apply the formula e A+B = e A e B e~^ A,B ^ which is valid 
for operators having c-number commutator. This gives 



2 
x-y 



E x k *k \ E * k+lx ^ E i E y k+l x>>,< 

r£^k — l g k,l — l — — 1 g k,l—l 



\ E * k ** I E * fe+! *<w / oo \ E ^ I E y k+l *k,i 

1 + E y m (A* + X m ) e*=i e w=* + 1 + E £ n (A> + X„) e*=* e fe -'= 1 (95) 

\ m=l / V n=l / 

The right hand side of this equation can be rewritten as follows 

oo oo 

1 + E my m [X m + E s n X n , m e- 1 «-.'=i + 

m=l \ n=l J J 

OO OO 

oo / _ oo ~ \ \ E y mx ™+2 E v k+ ' x kd 
l+J2nx n [X n + E y m Xn,m) e— 1 *.«=* . (96) 

71=1 V 771=1 / / 
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Due to (|35|) we obtain 



fc,i=0 ,3 fc,Z=0 



x—y 



\ 



( oo _ \ I E / oo _ \ I E 2/ fc+! ^,i 

E mx n i/ m X nim e k > l =° + E nx"y m X n , m e fc -'=° = 0, (97) 

yn,m=0 / \n,m=0 J 

what can be written in terms of the generating function fl3"7| ) as in (53). 
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